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Abstract 

An extension of the lower-bound lemma of Boggio is given for the weak forms of 
certain elliptic operators, which have partially Dirichlet and partially Neumann boundary 
conditions, and are in general nonlinear. Its consequences and those of an adapted Hardy 
inequality for the location of the bottom of the spectrum are explored in corollaries wherein 
a variety of assumptions are placed on the shape of the Dirichlet and Neumann boundaries. 
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I. Introduction 

E.B. Davies has promoted the use of operator and quadratic-form inequalities to 
obtain spectral information about Laplacians and Schrodinger operators. Among the goals 
of his analysis has been to understand the effect of the shape of the boundary on eigenvalues 
and the eigenfunctions. The present note may be viewed as a scholium to parts of the first 
chapter of [Dav89], where linear elliptic operators are bounded below, in the weak sense, 
by functions related to the shape of the boundary of a domain. The aim is the elaboration 
of lower bounds of similar kinds, given two complicating features: boundary conditions 
other than of Dirichlet type, and nonlinearities in the principal part of the operators. 

As has been pointed out in [FHT99], a precursor to many inequalities used for lower 
bounds in the spectral theory of elliptic differential equations is to be found a 1907 article by 
T. Boggio [Bog07] . This direct consequence of the Gaufi-Green theorem was interpreted 
in [FHT99] as a variational principle for the Laplacian and its nonlinear generalisation 
thep-Laplacian, on bounded domains with vanishing Dirichlet conditions on the boundary. 
Lemma 1.1 below, the main tool used in this article, likewise applies to a family of nonlinear 
elliptic operators, and non-Dirichlet boundary conditions are allowed. No great originality 
is claimed for Lemma 1.1: The idea is still close to that of [Bog07], and it is likely that an 
equivalent result could be extracted from what has been reported in [Maz81] or [OpKu90]. 
It is, however, stated in a form allowing geometric spectral results to be extracted easily, 
with attention to non-Dirichlet boundary portions. The short proof is given so that the 
exposition is self-contained. 

The great majority of lower bounds to the spectrum of elliptic partial differential equa- 
tions suppose Dirichlet conditions throughout any boundary that is present. One reason 
for this is that if Neumann conditions occur, the spectrum of the Laplacian is highly un- 
stable with respect to perturbations of the boundary. An arbitrarily small perturbation of 
the Neumann part of the boundary may cause the lowest eigenvalue to decrease arbitrarily 
close to 0. Intuitively, this is possible because the lowest eigenvalue of any domain with 
pure Neumann conditions is 0, so if a small such domain is weakly coupled to O via a thin 
neck or partitions, then the lowest eigenvalue will remain tiny. The instability of Neumann 
spectra has been extensively investigated, and instructive "rooms and passages" models 
date from Courant and Hilbert: See [CoHi37], [Maz73], [Fra79], [EvHa89],[HSS91],[Arr95], 
[BuDa02] among many treatments of the Neumann spectrum. 

There is evidence that if assumptions on the Neumann part of the boundary prevent 
its near-isolation, then the effect on the lowest eigenvalue will be more moderate. For in- 
stance, in [BuDa02], a continuous perturbation theory is established for purely Neumann 
Laplacians with respect to uniformly Holder perturbations. Section II of this article in- 
cludes lower bounds controlling the downward shift in the fundamental eigenvalue upon 
non-perturbative enlargement of the domain. 

It will be assumed throughout that O is a bounded domain in M d , the boundary of 
which is sufficiently regular that Green's formula is valid for suitable test functions and 
that the outward normal vector, denoted z/, exists a.e. at any part of the boundary where 
conditions other than vanishing Dirichlet are imposed. Regularity depends on the validity 
of certain Sobolev embeddings. Here each domain O will be assumed to have a boundary 
consisting of two pieces each of positive d — 1-dimensional Lebesgue measure, denoted 
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dVto and dVt^ := <9fi \ dQ^. For the purposes of this article it suffices to suppose that 
the boundary is of "class C" (see [EdEv87], p. 244, 248), and that the boundary is C 1 
a.e. with respect to the surface Lebesgue measure. Such a domain will be termed regular. 
(Regularity is not actually required on dVto, while on dVt^ usefully wider conditions of 
regularity may possibly be gleaned from [Ken94] . The widest conditions for the boundary 
shall not be pursued here.) An adapted set of test functions V (fi, £Id) can be defined as 

consisting of the restrictions to O of 

Although the notation <9fi n is introduced with Neumann boundary conditions in mind, 
no such assumption is made. Only the Dirichlet portion of the boundary is actually subject 
to (weakly defined) boundary conditions. This article is concerned solely with functionals 
on V, without investigating the domains of definition of the elliptic operators that the 
functionals may define. 

Lemma LI. Suppose that O is regular, 1 < p < oo, Q is a real-valued absolutely 
continuous vector Geld, and A is a real-valued absolutely continuous tensor Geld of type 
Ti . (In elementary terms, A is a d x d real matrix-valued function acting on Q.) Suppose 
further that v ■ A T Q < a.e. on Q N . Then for all ( e V (Q,Q D ), 

JjAVCfdx > J jdiv (A T Q) - (p-l)|Q| p '}|Cr dx, (1.1) 

where p' is the dual index, p' = . 

Alternatively, Let a and q be absolutely continuous functions on O and let u be a fixed unit 
vector such that aqu • v < a.e. on 8Qn- Then for all £ EV (O, Qd), 

J \au-V(\ p dx> J {u-V(aq)-(p-l)q p '}\(\ p dx. (1.2) 

Proof. We calculate following [FHaT99] with the divergence theorem: 

< - / v ■ (A T Q|Cr) dS = - f div (A T Q|cr) dx 

JdQ JO. 

= - [ (divA T Q)\C\ p dx-p [ \(\ p - 2 (Q-AV(dx. 
Jn Jet 

With w = |C| P_2 CQ 7 the arithmetic-geometric mean inequality [HLP34] gives 

i (Avo ■ w| < - \awc\ p + 1 |wK = - \awc\ p + ^ icriQK- 

p p p p 

Collecting terms, 

J \AV(\ p dx > J (div (A T Q) - (p - 1) \Q\ P ') \(\ p dx, 
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which is (1.1). The proof of (1.2) is strictly analogous, beginning with the Gaufi-Green 
formula in the version 



< - 



(aq\C\n SdS = - 



an 



d(aq\(\P) 



an 



dxi 



dx 



for a Euclidean coordinate Xi (e.g., [Ev98]). The result is 



* a, 



a 



dxi 



dx > 



\C\ p dx, 



which when written without coordinates becomes (1.2). □ 
Remarks 

1. Inequality (1.1) is well-known in the case p = 2 with Dirichlet boundary conditions 
on the entire boundary, in the form obtained with 

Q std := -AVlog($), 
for a suitable function $ > on O. For p ^ 2 the analogous choice is 



Q 



std 



|AV$| P " 2 AV$ 



A simple calculation shows that for this choice Lemma 1.1 becomes 

-div (|AV$r~ 2 A T AV^j 



\AV(\ p dx > 



n 



n 



| C |i 



dx, 



(1.3) 



which for p = 2 is essentially Theorem 1.5.12 of [Dav89]. Let us refer to this as the 
standard form of Inequality (1.1). It is all the more familiar for the linear Dirichlet 
Laplacian, i.e., with A taken as the identity and p = 2, in which case it becomes the 
well-known bound 

which is often attributed to Barta or Duffin, who, however, published much later than 
Boggio [Bog07]. (In the not widely accessible [Bog07], Boggio used the Gaufi-Green 
Theorem to prove the case of Lemma 1.1 with d = 2 = p, pure Dirichlet, A = the 
identity.) Lemma 1.1 is truly more general than the standard form, which corresponds 
to a vector field Q restricted to be A times something irrotational. 
2. Similarly, if 



Qstd 



' dxi 



p-2 



dxi 



-1 



is inserted into (1.2), the result is 



n 



,3C 

dxi 



dx > 



d 

dxi 



' dxi 



p-2 



dxi 



-l 



\(\ p dx, 



(1.4) 
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provided that aqv 1 < a.e. on dVL N . The expression in the numerator of the lower 
bound (1.4) is a one-dimensional p-Laplacian. 
3. Lemma 1.1 may be thought of as a lower-bound variational principle for the bottom 
of the spectrum of the operators whose energy forms are as in the left side of (1.1), 
i.e., 



If $ > is the fundamental eigenfunction of H, that is, if$ = Ai$ p_1 , then the lower 
bound to the spectrum of H thus obtained reduces to Ai. In this sense the inequality 
is optimal. While it must be conceded that the inequality is thus arguably not better 
in theory than the standard form for the purpose of finding lower bounds, it has two 
practical advantages: 

a) It is relatively easy to construct a vector field with positive divergence to cover a 
domain the only material restriction being that it be incoming on 8Qn. 

b) The nonlinearity of the bound can be exploited. If any vector field can be con- 
structed on Q with the necessary boundary behavior and strictly positive diver- 
gence, then by scaling Q — > tQ a strictly positive lower bound will be obtained 
for sufficiently small t > 0. 

4. An extension from Euclidean domains to subdomains of orientable smooth manifolds 
is straightforward. 

In the case where A is the identity tensor, H corresponds to the p-Laplacian, denoted 
— A v . It is known that the minimum of the functional 



for \\C\\p — 1 and Dirichlet boundary conditions, corresponds to an eigenvalue Ai of — A p 
Many other familiar facts for the linear case p = 2 carry over, such as that the minimiser 
u\ exists in the appropriate Sobolev space and is positive on O. The second eigenvalue of 
(1.5) is also inflection point of (1.5) and can thus be characterised variationally. However, 
the theory of the higher eigenvalues of — A p remains murky when p ^ 2. Detailed analysis 
of the p-Laplacian will not be entered into here, as the focus will be entirely on lower 
bounds for its energy form H. (The spectral theory of linear elliptic differential operators 
is presented, for example, in [EdEv87], [Dav95], and a useful review of the spectral theory 
of the p-Laplacian is to be found in [DKN96, DKT99].) 

In the following section particular vector fields Q will be chosen responding to as- 
sumptions about the shape of the domain O. 

II. Spectral bounds with Dirichlet conditions on subsets of the boundary 

The aim of this section is to find lower bounds to the spectrum in terms of the shapes 
of 8Qd and dQjy. For simplicity it will be assumed henceforth that A is the identity tensor, 
i.e., the case of the p-Laplacian. No material difficulty would arise if A were retained. 

In essence this section consists of a selection of examples of the use of Lemma 1.1, 
organized as a list of corollaries. 




(1.5) 
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Definition II. 1. Let Ai (Q,dQr>) denote the infimum for £ e T> (Q, 8Qd) (and not iden- 
tically zero) of 

'Liver 



E P (0 



In icr 



The quantity Ai (fi, dVtc) will be referred to as the fundamental eigenvalue of the p- 
Laplacian with respect to fi, dVto- 

As remarked in the introduction, eigenvalues can be rather unstable with respect to 
Neumann conditions on the boundary. In the purely Dirichlet case the principle of domain- 
monotonicity is familiar: If fi C T, then Ai (fi,pureDBC) > Ai (T,pureDBC). There is 
no such general principle in the presence of Neumann conditions. Nonetheless, Lemma 1.1 
implies the following: 

Corollary II. 2 — Restricted Monotonicity Principle. Let fi and T be regular do- 
mains such that fi C Y. Let ui(T,8Td,p) > on fi be a fundamental eigenfunction 
for the p-Laplacian with respect to T, T^. If v ■ Viti(T, dYp^p) > a.e. on 8Qn, then 

\ 1 (n,dn D ,p)>x 1 {r,dT D ,p). 

Proof. The assumptions allow the standard form (1.3) of the lower bound with the choice 
$ = iti(T, dT D ,p). In this case the right side of (1.3) reduces to Ai (T,dT D ,p) . □ 

An immediate consequence is that if the Neumann boundary is a graph with respect 
to a hyperplane, there is an analogue of the elementary outradius bound for the Dirichlet 
problem. There results a crude "box bound" in terms of fii :P : = the fundamental eigenvalue 
of the one-dimensional p-Laplacian — Ap on the unit interval, with Dirichlet conditions 

2 

at one end and Neumann at the other. (Thus = i - -) 

Corollary II. 3 — Box Bound. Let Q be a regular domain such that dQjy is the graph of a 
function on a subset of {x : Xi = 0}, with i/-e, < a.e. on dQ n For L := sup (xi : x G fi) — 
inf (ii : x 6 Q), 

Ai (Q,dQ D ,p) > 



fJ>i, P 



LP 

Proof. Since the fundamental eigenfunction ui jP (x) of the one-dimensional p-Laplacian 
with Dirichlet conditions at one end and Neumann at the other is positive and monotonic 
(due to the maximum principle), we may choose a = 1, $ = ui jP in (1.4), yielding the 
result. □ 

Less crude bounds involving the shape of O can be obtained with more sophisticated 
comparisons. Recall that in [Dav98], Davies exploited Hardy's one-dimensional inequality 
[Har20] [HLP34] to produce lower bounds in the quadratic-form sense of the type 

-A> ° 



(m(x)) 2 ' 



where m is an averaged distance to the boundary, supposed entirely Dirichlet. These 
bounds extend almost immediately to the case of the p-Laplacian and boundaries that are 
only partly Dirichlet, with certain restrictions: 
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Definition II. 4. Given a regular domain Q, the Dirichlet boundary sector of x, dfi^, 
denoted S 1 (x, dQr>), is the union of all line segments contained in O joining x to dQ,o- 
Following [Dav89], an averaged distance to the Dirichlet boundary is given by 

1 1 r dS(u) 

(m(x)f : " ~~ d J ' K (rf u (x)f 

where the measure dS(u) is the Lebesgue measure on the unit sphere S ,d_1 , ua is the total 
measure of K := S (x, dflo) H and d u is the length of a line segment joining 

x to dflo- 

Note: The Dirichlet boundary sector of x may be vacuous, and ^ may be zero. 
Lemma II. 5. Let O be a regular domain and £ e X> (f2, D ). Then for any 1 < p < oo, 



/ |VCrda;>p- p (p-l) p - 2 / 

JO, J Q 



1 


CM 


Iq. 


m(x) 



With the one-dimensional L p Hardy inequality (i.e., [OpKu90] Theorem 1.14 with p = q, 
v = 1, and w(x) = (x — a) _p , itself derivable from lemma 1.1), the proof of the lemma is 
almost exactly as for [Dav98], Theorem 1.5.3, and will be omitted. 

Suppose that <9Qat is starlike with respect to an exterior point, at which the origin will 
be placed. It will now be shown that there is a sort of multidimensional Hardy inequality 
with respect to a distance function that is the intuitively the lesser of a multiple of |x| and 
the distance to the Dirichlet part of the boundary. The precise version of this statement is 
expressed by (2.2) in Corollary II. 6 and (2.3) in Corollary II. 7. Versions of these corollaries, 
in the linear case p = 2, have been previously reported, though not published, by the author 
since [Har93]. They imply, for example, that it is possible to have a fractal 8Qn without 
the collapse of the bottom of the spectrum to 0. 

Corollary II. 6. Let 1 < p < d, and suppose that the origin is exterior to O, that 8Qn 
is starlike with respect to the origin, and that the interior of O adjacent to dfl n is away 
from the origin. Then for all C G V (O, 90^, p), 



|VC| P ete> ( d -^' 

Q V V 



C 



dx. 



(2.1) 



Consequently, 

Ai {n,dn D , P ) > 

Moreover, for any 7, < 7 < 1, 



d — p 



pinf {|x| : x ^ 0} 



A 1 (0. iKln.p) > inf I 7 ( Vrf + & ~ 

p-i 



(2.2) 



Proof. Eq. (2.1) results from the choice Q = (y^rj f^p" m Lemma 1.1 and a cal- 
culation. The next statement is the immediate lower bound by replacing a factor in the 



7 



integrand by its infimum. Statement (2.2) is obtained in the same way from the weighted 
average of the lower bounds of Lemma II. 5 and (2.1). □ 

When p > d, the origin needs to be taken from the other side of 8Qn- 

Corollary II. 7. Let d < p < oo, and suppose that the origin is exterior to O, that 8Qn 
is starlike with respect to the origin, and that the exterior of adjacent to d£l n is away 
from the origin. Then for all £ E V (fi, dVlrj,p), 



L mPdx -{ 



Consequently, 

X 1 (Q,dQ D ,p) > 
Moreover, for any 7, < 7 < 1, 



Ai (fi, dVtu.p) > inf 
xef2 
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p — d 
p 



c 



dx. 



(2.3) 



V 



d 



p'vai {|x| : x ^ 0} 



p — d 
p|x| 



+ (1-7) 



(p_ l)P-5 

p p m('x) p 



(2.4) 



Proof. Eq. (2.3) results from the choice Q = - 
calculation. The other statements follow as in Corollary II.6. 



p-i 



in Lemma I.l and a 
□ 



For the missing case p = d, a comparison can be made with annular domains, i.e., the 
region between two concentric circles or spheres. 

Definition II. 8. Let a(r,R,p,d) denote the lowest eigenvalue of the p-Laplacian on 
A := {x : r < |x) < R}, with Dirichlet boundary conditions imposed where |x| = R 
and Neumann boundary conditions imposed where |x| = r. Similarly, let /3(r, R,p,d) 
denote the lowest eigenvalue of the p-Laplacian on A, with Neumann boundary conditions 
imposed where |x| = R and Dirichlet boundary conditions imposed where |x| = r. 

Observe that ct(r, R,p,d) and f3(r, R,p, d) can be written in terms of explicit Bessel 
functions, and that for any p > 1 a and (3 are determined by ordinary differential equations, 
as the fundamental eigenfunctions are radial. The eigenvalues ct(r, R, 2) and f3(r, R, 2) and 
their associated eigenfunctions are thus numerically accessible. The Restricted Monotonic- 
ity Principle II. 2 immediately implies: 

Corollary II. 9 — Annulus Bound. 

a) Let p > 1 and suppose that the origin is exterior to O, that dQpj is starlike with 
respect to the origin, and that the interior of O adjacent to 8Qn is away from the 
origin. Suppose further that O C {x : r < |x| < R}. Then Ai (Q, dVt,p) > a(r, R, p, d). 

b) Let p > 1 and suppose that the origin is exterior to O, that 90^ is starlike with 
respect to the origin, and that the exterior of O adjacent to 90 n is away from the 
origin. Suppose further that O C {x : r < |x| < R}. Then Ai (O, dfl,p) > f3(r, R,p, d). 
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Remark. The Annulus Bound, as well as the Box Bound II. 3, can be combined with 
Lemma II. 5 to produce lower bounds analogous to (2.2) and (2.4), taking the Dirichlet 
boundary into account. 

Finally, let us briefly consider the effect of changing a portion of the boundary from 
C <9fijv to C 8Qd or vice versa. To focus the analysis, it will be supposed that O is 
convex, and that the unperturbed problem has boundary of one type, either <90 = <90tv 
or <90 — dO,o- 

The former case is easier to analyze, and indeed the perturbation theory when ad- 
ditional Dirichlet conditions on small subsets of positive capacity has been studied sys- 
tematically, for example in [Oza81], [Oza82], [Cou95], [McG96], [McG98]. A rough non- 
asymptotic estimate follows from Lemma II. 5: 

Corollary 11.10. Let d > 2, 1 < p < oo, and suppose that O is strictly convex. Define 
e(0, p) as the minimum for x, y G <90 and v y a unit normal to a support plane at y of 
|yZx|p+'* ' (This is a measure of the eccentricity of O in relation to its diameter.) Then 

Ai (o, ao D , P , d) > P -v{p - 

Proof. Strict convexity implies that a change of variable can be made in the integral 
in Definition II. 4 of -, — tU-w passing from u G K C S d ~ l to y G 90^, as K is the 

projection to the unit sphere of 8Qd- This requires a Jacobian factor ^ , v y being 

|y x ! 

the outward normal where defined. (By assumption in this article the outward normal is 
defined a.e., but it is convenient to let u y range over the normals to support planes, for 
definiteness at all boundary points and so the formula will remain valid for some irregular 
domains.) The infimum of the transformed integrand is decreased by allowing x and y to 
vary independently over the closure of O. In that circumstance, it is easy to see that the 
minimum is attained for x, y G 0, yielding the result. □ 
The problem of changing a portion of a Dirichlet boundary to Neumann is more del- 
icate than the converse. If a portion of the boundary is changed from having Dirichlet 
conditions to Neumann, then a positive lower bound to Ai can be exhibited under starlike- 
ness conditions of Corollaries II. 6, II. 7, or II. 9. Such bounds will, however, fail to approach 
the unperturbed eigenvalue as |<90jv| — > 0. The perturbation theory for the introduction 
of Neumann boundaries has apparently been studied only under very special conditions 
to date (e.g., [Oza85]). A full theory will probably require a careful consideration of the 
asymptotic behavior of eigenfunctions near the perturbation. 
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